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We calculate entanglement entropy in a non-relativistic field theory described by the Schrodinger 
operator. We demonstrate that the entropy is characterized by i) the area law and ii) UV divergences 
that are identical to those in the relativistic field theory. These observations are further supported 
by a holographic consideration. We use the non-relativistic symmetry and completely specify entan- 
glement entropy in large class of non-relativistic theories described by the field operators polynomial 
in derivatives. The entropy of interacting fields is analyzed in some detail. We suggest that the 
area law of the entropy can be tested in experiments with condensed matter systems such as liquid 
helium. 
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Introduction. Entanglement entropy of a field system 
is defined by tracing over modes that reside inside a 
chosen surface E If the total system is character- 
ized by a pure quantum state the subsystem inside the 
surface is described by a density matrix p with the von 
Neumann entropy S = — Trplnp known as entangle- 
ment entropy (see recent reviews Q , [I| and references 
therein) . This entropy is non- vanishing provided there 
are short-distance correlations in the system. The 
presence of such correlations has two consequences: 
a) entanglement entropy is determined by geometry 
of the surface E, to leading order by the area of E and 
in higher orders by intrinsic and extrinsic geometry of 
E, see and b) the entropy is UV divergent. In d 
space-time dimensions, for a field system described by 
a relativistic wave equation, one finds that 

s~iv^P. (i) 

The exact coefficient of proportionality in (JTJ depends 
on the regularization scheme, N is the number of held 
species. Entanglement entropy has been actively stud- 
ied in the literature. The main focus has been made on 
the relativistic field systems and especially the confor- 
mal field theories. In the latter case a holographic de- 
scription of entanglement entropy has been proposed 
that suggests an alternative way to compute the en- 
tropy in the strongly coupled conformal field theory 
dual to supergravity on anti-de Sitter space-time Q . 

The definition of entanglement entropy is applicable 
to any system that has a quantum mechanical descrip- 
tion and does not necessarily tie to the Lorentz sym- 
metry. A typical example considered in the literature 
is a finite system of coupled non-relativistic oscilla- 
tors [![. The Hamiltonian of this system is quadratic 
in momenta. The continuous limit then reproduces a 
relativistic field described by the Klein-Gordon opera- 
tor which is quadratic in time derivative. The primary 
purpose of the present note is to compute the entan- 
glement entropy for a field system described by the 
Schrodinger field operator linear in time derivative. 



The field system in this case can not be presented as 
a collection of harmonic oscillators. Nevertheless, as 
we show below, the structure of the entropy in this 
case is essentially identical to the one in a relativistic 
field theory. 

The AdS/CFT correspondence can be generalized 
to non-relativistic conformal field theories i (for an 
earlier work see @), see also Q- These theories de- 
scribe the "fermions at unitarity" . In a wider context 
the non-relativistic theories effectively describe many 
realistic systems well studied experimentally and the- 
oretically in condensed matter, liquid helium He 4 is 
the standard example. Our results thus are applica- 
ble to these systems too. 

The field models with the broken Lorentz symmetry 
are sometimes considered as UV completions of other- 
wise non-renormalizable relativistic effective theories, 
e.g. the General Relativity. We anticipate that our 
results can be extended to those theories. 

The replica method. Before proceeding we remind the 
technical method very useful for calculation of en- 
tanglement entropy. This method is known as the 
replica method, see ref.Q. One first observes that 
— Tr phi p = (ad a — l)Tr p a \ a =i- The next observation 
is that the density matrix obtained by tracing over 
modes inside the surface E is Trp Q = W[a], where 
W[a] = — hxZ[a) and Z(a) is the partition function 
of the field system in question considered on Euclidean 
space with a conical singularity at the surface E. Thus 
one has that 

S=(ad a -l)W(a)\ a=1 . (2) 

One chooses the local coordinate system {X^ = 
(t, afj)}, where r is Euclidean time, such that the sur- 
face E is defined by conditions r = 0, Xi — and 
(x2, ■■,Xd) are the coordinates on E. In the subspace 
(t, xi) it is convenient to choose the polar coordinate 
system t = r sin(0) and x\ = r cos((/>) where angular 
coordinate <j> changes in the limits < <j> < 2tt. The 



2 



conical space in question is then denned by making 
the coordinate <f> periodic with the period 27ra, where 
(1 — a) is very small. 

In order to calculate the effective action W(a) we 
use the heat kernel method. For manifolds with coni- 
cal singularities this method was earlier developed in 
great detail in [9J. Consider a quantum bosonic field 

i/2 

described by a field operator T> so that Z = det ' V. 
Then the effective action is defined as 



W = - 1 - r*L<*K 

2 7,2 s 



where e is an UV cut-off, is expressed by means of the 
trace of the heat kernel K(X,X',s) = (X\ e - sV \X') 
satisfying the heat kernel equation 



(d s +V)K(X,X',s) = , 
K(X,X',s = 0) = 5(X,X') 



(4) 



In the Lorentz invariant case the heat kernel 
K(cf>, <f)', s) (where we skip the coordinates other than 
angle <f>) on regular flat space depends on the differ- 
ence {(f) — <f>'). The heat kernel K a ((f>,<f>,s) on space 
with a conical singularity is then constructed from this 
quantity by applying the Sommerfeld formula 10] 

K a (<f),<t>', S )=K{cf>-<f)',s) 

% f w 

+- / cot — K{6 - <f>' + w, s)dw . (5) 

Ana J r 2a 

The contour T consists of two vertical lines, going from 
(— n+ioo) to (— 7T— too) and from (it— too) to (n — Hoo) 
and intersecting the real axis between the poles of the 
cot —2ira, and 0, +2ira respectively. For a = 1 
the integrand in (J5J) is a 27r-periodic function and the 
contributions of these two vertical lines cancel each 
other. Thus, for a small angle deficit the contribution 
of the integral in ([3]) is proportional to (1 — a). 

The Laplace operator. In d space-time dimensions, for 
a relativistic theory described by the massless Klein- 
Gordon operator or, in Euclidean signature, by the 
Laplace operator V 2 = d 2 + X)j=i ®1i where r is Eu- 
clidean time, the heat kernel is known exactly 



K(t, t , x, x , s) 



1 



( 47rs )d/2 



(6) 



We take (d — 2)-surface E to be the plane defined 
by equations x\ = 0, r = so that {x%, £3, .., x^) 
are coordinates on E. In the polar coordinate sys- 
tem t — r sin <f> and x\ — r cos <f> we have for two 
points (r, cf>) and (r, </>') that (r - t') 2 + (x\ — x[) 2 = 
4r 2 sin 2 (^^). The trace is defined as Tr K a = 

J d d ~ 2 Xi J °° dr r J^ a d4>K a (4> = cf>',r' = r^x t = 
Xi,s). For the contour integral T we find that 

sy , \ i f ^ w dw 1 2 

C 2 (a) = — / cot — = — 1 - a 2 . 7 

8na J r 2a sm ^ 6a 2 



Thus one obtains for the trace of the heat kernel (j^j) 
1 



TvK a = 



(4irs) d /' 



{aV + s 27raC 2 {a)A(^)) , (8) 



where V = J drd d 1 x is the volume of spacetime and 
A(E) = J d d ~ 2 x is the area of surface E, and 



2 tt ACE) 

{d) 3(d - 2) (47r) d / 2 e d - 2 



(9) 



(3) for entanglement entropy in d space-time dimensions 



The Euclidean Schrodinger operator. We consider a 
non-interacting system described by a non-relativistic 
field operator V — —2imd t — Y^tZid 2 - Written in 
this form the Schrodinger operator T> can be viewed 
as a reduction of the massless Klein-Gordon operator 
in d+1 dimensions 2d^dt^df after projecting onto the 
space of fixed momentum — —im as described in 
0, @. The Euclidean Schrodinger positive operator 
is further defined as 



V E = 2m y ^d 2 - df 



(10) 



which is formally obtained by a double analytic con- 
tinuation t — y zt, m — y —im. The heat kernel of this 
operator is a product of the heat kernels of two oper- 
ators 

s - s (- 9 ?) 

The heat kernel of the (d — l)-dimensional Laplace 
operator &f is given by expression similar to ([B]). As 
for the heat kernel of the first order operator 2vmd T it 
is convenient first to Laplace transform [TTJ] it to the 
heat kernel of the second order operator — 9 2 , 



-s2m\ 



da 



III 



vW/ 2 



(11) 



and then use expression of the type ® for the heat 
kernel of operator — d 2 .. Switching to the polar co- 
ordinates (r,<p), t — rsin0, x\ — rcos^, we find for 
the heat kernel function of two points characterized by 
different angular coodinates <f> and <j)' (r' ~r,x[ — Xi), 



b'le'^ 



E U\ — 



2ms 



(47TS)~ 

f(^™s 2 ± + l S in 2 ±) 



where we introduced x 



and if) 



2 /° 



(12) 



The kernel now depends on two angular variables 
if) and x- This is a consequence of the absence of the 
Lorentz symmetry. Constructing the heat kernel on 
space with a conical singularity we have to make the 
heat kernel to be periodic function of both angles, x 
and if), with the period 2ira. We are however inter- 
ested in the result to leading order in (1 — a). The 
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Sommerfeld formula applied to the angular coordinate 
X already gives the term proportional to (1— a). Thus, 
it is sufficient to apply the formula (O only to the 
angular coordinate \ an d consider the other angular 
variable tj) as periodic with the period 2tt. Calculat- 
ing the trace one thus first replaces \ by X + w , takes 
4>' = 4>, integrates over <p from to 27r and then applies 
formula |5| to w. Using the integral 



which is invariant under the scale transformations 



2tt 



do 



drre 



f (icos 2 0+i; 



sin 2 f 



and J °° d<r<j- v e-( sm f/ a = (sm) 2 ^-^T{v - 1) and 
formula (0 we arrive at the expression for the trace 
of the heat kernel 



Tr e 



-sVb 



2m~ 1 aV 

(47TS) 2 



(47TS) 



T 2^0 2 (a)A(£) , (13) 



valid up to terms of order (1 — a) 2 . We see that the 
volume part of the trace of the heat kernel has the de- 
pendence on the proper time s typical for the heat ker- 
nel of the Laplace operator in dimensions. This 
seems to be consistent with the way we obtained the 
Schrodinger operator as a reduction from the Klein- 
Gordon operator in (d + 1) dimensions. 

On the other hand, rather surprisingly, the surface 
term in (|13| is precisely identical to the surface term in 
the heat kernel of the d-dimensional Laplace operator 
p|). The entanglement entropy in a non-relativistic 
theory in any dimension d is thus given by the same 
formula © as in the case of a real scalar field de- 
scribed by the Lorentz invariant Laplace operator. In 
four space-time dimensions we hence obtain 



5, 



(d=4) 



487T6 2 



(14) 



In usual signature the non-relativistic field is natu- 
rally described by a complex function. For a complex 
field everything we have just calculated then should 
be multiplied by 2. 

The fact that the entropy does not depend on the 
mass parameter m has a simple explanation. The op- 
erator T>e is invariant under the rescaling t — > /3t, 
m — > (3 m. Note that both terms in f|13|) are invariant 
under this rescaling. Since the entropy is not supposed 
to be a function of r (or t in the usual signature) it is 
not a function of m either. 



A holographic consideration. Recently there have 
been proposed some dual gravity descriptions of non- 
relativistic systems. In this note we use the ap- 
proach suggested in Q. As an example we consider 
5-dimensional geometry described by metric 



dr 2 

ds 2 = L 2 (^ 



1 



d^_. 



(15) 



— > Xxi 



Xr 



(16) 



For z — 2 these are also the symmetries of the 
Schrodinger operator T>. Thus, one expects that (for 
z = 2) the metric (|15[) gives us a gravity description 
of the non-relativistic system. For z = 1 metric (fT5"j) 
is the anti-dc Sitter metric which has larger symmetry 
50(4,2). In this case the AdS/CFT correspondence 
is at work with a precise dictionary of translation be- 
tween the quantum field theory and the gravity de- 
scriptions. 

In the AdS/CFT correspondence it is known that 
entanglement entropy has a holographic gravity de- 
scription (for a recent review see [3j). The entropy is 
defined in the quantum theory living on the boundary 
of the 5-dimensional space-time (fl5|) . Since the geo- 
metric quantities are usually divergent when extended 
till the boundary one considers a regularized bound- 
ary at r = e. According to the AdS/CFT dictionary, 
e plays the role of the UV regulator on the field the- 
ory side. For a given 2d-surface S on the regularized 
boundary consider a minimal 3d-surface 7 such that 
it lies entirely on the slice t = const and dj — Y,. 
According to the holographic proposal the quantity 



Shol — 



Areafa) 
4G 5 



(17) 



computed on the gravity side with Newton's constant 
G5 , is identical to entanglement entropy computed for 
the quantum fields defined on the 4d-boundary of the 
space-time (jT5j • The correspondence is easily estab- 
lished for the UV divergent part of the entropy as soon 
as we note that 



Area{i) = ^ A &) 



(18) 



where .. stands for the subleading terms. Moreover, 
in the AdS/CFT correspondence one has a relation 



TIL 3 
205 



= N 



(19) 



where N is the number of fields (species) in the bound- 
ary theory. In fact, the theory on the boundary is a 
SU(n) Yang-Mills theory, n = y/N is the number of 
colors. Thus, using the holographic proposal one ob- 
tains 



S h oi = N 



MX) 

47T6 2 



(20) 



i=l 



which is the right expression for entanglement entropy 
of N relativistic fields in four space-time dimensions. 
(Note that the holographic UV regulator e is not iden- 
tical to the proper time regulator e appearing in the 
heat kernel calculation.) 

We now notice that the induced metric on a slice of 
constant t of ([T5]l does not depend on the parameter z. 
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Thus, for a given surface E on the regularized bound- 
ary the minimal surface 7 and its area are the same 
for any value of z. Assuming that the holographic de- 
scription of entanglement entropy can be extended to 
other values of z we obtain that the entropy then is 
still proportional to the area of E and is UV divergent 
in the same way as for the relativistic fields. For z = 2 
the theory on the boundary is a non-relativistic the- 
ory described by the Schrodinger operator. The holo- 
graphic description thus is consistent with our direct 
calculation of entanglement entropy in a field theory 
with the Schrodinger operator. 

Furthermore, from the fact that neither of the quan- 
tities that appear in (|17|) depends on z it is natural to 
conclude that relation (TK))) still holds for value z = 2 
so that N in this case is the number of degrees of free- 
dom in the non-relativistic theory. The holographic 
formula (T30|) then is consistent with the findings in 
this note: N non-relativistic fields produce same en- 
tanglement entropy as N relativistic fields. Clearly, 
this statement and the equation (|2T)]) . can be extended 
to any dimension d. We note that in higher dimen- 
sion d entanglement entropy is given by expression 
([9]) (the exact numerical prefactor may depend on the 
regularization scheme) and is invariant under the scale 
transformations ([T5]) . For z = § the holographic en- 
tropy was earlier calculated in [l2j and agrees with 
our result. 

We also note that the invariance under transforma- 
tions (fTo]) with A = 2 helps to understand the struc- 
ture of the heat kernel (fl3|) . Under (fTB]) the proper 
time s rescales in the same manner as r 2 , s — > X 2 s. 
Since the volume V changes as V — > X d+1 V one needs 
s~5~ in the denominator of the volume term in (|13l) 
for the invariance under transformations (TTJ 



entropy of interacting fields in d dimensions 

'(<*) 



*../, = *^ a Ap)M*-T(l - %M 2 e 2 ) , (21) 
12(47r)— *> 



where we used that s d l 2 e Aj2s 



ra- 



|,MV). 

Clearly, the leading UV divergence of the entropy (12T1) 
is again (multiplied by N) (j9|) and is thus not affected 
by the presence of the interaction in the action. This 
is consistent with the holographic calculation. As we 
have seen, the holography predicts that the leading 
UV term in the entropy (|20|) is the same as for free 
fields. 

The interaction however shows up in the sub- 
leading UV divergent and UV finite terms. For in- 
stance, in four dimensions we find 



S = 



NA(H) 1 



48tt 



-2 + M% c )( 7 - 1 + ln(eM(^ c ))^)](22) 



valid both in relativistic and non-relativistic cases. In 
the holographic description the subleading terms in 
the entropy appear as the subleading terms in the area 
of minimal surface (|f 8p . The parameter M in (|2"21 
depends on the (constant) background field ip c so that 
one may apply the renormalization technique similar 
to one developed in [l3[ for cosmic strings. We also 
notice the earlier work on entanglement entropy in 
relativistic O(N) model [l~4|. 



Polynomial field operators, symmetry and the entropy. 
The symmetry argument helps to obtain the general 
structure of the entropy in the case when the non- 
relativistic theory is described by a polynomial field 
operator 



-d 2 T + m 2(1 - n H-d?) n 



(23) 



The heat operator (d s + T>p) in this case is invariant 
under the transformations 



The interacting fields. So far we have considered the 
non-relativistic fields without interaction. The holog- 
raphy is however supposed to describe a strongly cou- 
pled field system. In order to check the holographic 
predictions we have to introduce interaction. The in- 
teraction can be included by adding a potential term 
j d d xV(ip) to the classical action, here <p is a set of 
fields in question. 

In the one-loop approximation one splits <p — ip c + 
ip q , where ip c is the classical background field and ip q 
is the quantum field. The integration over tp q then 
reduces to calculation of the functional determinant of 
operator V E + M 2 (ip c ), where M 2 = V"((p c ). If there 
are iV fields, for simplicity we assume that M 2 is the 
same for all fields. A generalization of this simplified 
situation is however straightforward. The heat kernel 
of this operator is the product e~ sT>E ■ e~ sM . Using 
the already calculated trace flU)) we obtain for the 



X 2n s 



and Xi — > (3xi , m — > /3 1 -™ m . (24) 
The trace of the heat kernel then should have the form 



Tr a e 



-sV P _ 



(d-lVn-l) 



S 2™ 



tV 



+bC 2 (a)(- 



,1/2 



.)^(E) , 



(25) 



which is invariant under (|24[) . where a and b are some 
constants which may depend on n and d. Entangle- 
ment entropy in a quantum field theory described by 
the polynomial field operator (|23|) is 



-)^A(E) 



(26) 



For n > I the UV divergence of entanglement entropy 
becomes milder than in the relativistic case ([T]). 
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Conclusions. Entanglement entropy is determined 
by the short-distance correlations in the field system. 
That's why it is surprising that the non-relativistic 
theory is characterized by the same entropy as a rela- 
tivistic field theory although the short-distance behav- 
ior of two theories is different. (Indeed, in d dimen- 
sions the equal time correlation function is l/|x| d_2 
for a relativistic scalar field and <5(x) for a non- 
relativistic field.) In the present note we have demon- 
strated this by a direct calculation of entanglement 
entropy using the replica method and by extending 
the holographic description of entanglement entropy 
to the non-relativistic theories which have dual gravity 
descriptions. The agreement between two methods, in 
particular, indicates that the holographic description 
of entanglement entropy should have wider applica- 
tions. 



The fact that we have the area law in the non- 
relativistic case may have some interesting conse- 
quences. In particular, the area law could be checked 
experimentally One would have to prepare a pure 
state in a condensed matter system available at the 
laboratory and then arrange a situation when part of 
the system is "hidden" for observations. This may be 
yet another way to find in condensed matter the fea- 
tures typical for black holes, the analog of the Hawk- 
ing radiation is the well known other example. 

This paper is based on the talk which the au- 
thor gave at "Symetries non relativistes : theorie 
mathematique et applications physiques" , Tours, 23- 
24 juin 2009. The helpful discussions with A. Barvin- 
sky, D. Fursaev, O. Lysovyi are kindly acknowledged. 
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